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SUHMAKY 


The  use  of  high  altitude  balloons  as  stable  platforms  for 
meteorological  and  astronomical  observatories  necessitates  ai  investiga¬ 
tion  or  the  dynamic  oenavior  ol  balloons.-  This  report  deals  with  tne 
azimuthal  rotations  of  the  balloon-gondola  system..  The  aerodynamic 
damping  and  inertia  of  the  system  have  been  mattiemat ical Iy  formulatea 
ana  experimentally  verified.  The  forcing  torque,  mostly  aerodynamic  in 
nature,  has  been  described  and  its  order  of  magnitude  has  oeer.  estimated 
for  known  balloon  rotations.; 


ft' 


These  rotations  have  been  alleviated  in  the  past  by  control 
systems  which  orient  the  gondola  by  reaction  torques  applied  to  tne 
balloon.-  The  aerodynamic  parameters  which  are  evaluated  m  tnis  report 
can  ba  useful  in  the  design  of  such  a  control  eystem. 


I*.  ’✓-! 
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II ■  INTRODUCTION 


We  believe  the'  the  principal  source  of  the  azimuthal  motion 
are  torques  generated  in  the  wake  of  the  balloon.,  as  noted  in  our 
Report  No.  1  of  this  contract  (Ref.  1).  As  the  balloon  moves  through 
the  atmosphere,  the  flow  of  air  past  the  balloon  separates,  forming  a 
cylindr. cal  vortex  sheet  in  the  wake.-  The  sheet  is  hignly  unstable  and 
curls  upon  itself,  thus,  creating  a  region  of  highly  concentrated  vor- 
ticity.  Eventually,  this  region  breaks  awey  from  the  sheet,  and  the 
process  is  repeated.  The  overall  result  is  a  periodic  shedding  of  vor¬ 
tices  in  the  wake.  Aerodynamic  torques,  due  to  these  vortices,  can 
Induce  the  balloon  to  exhibit  many  modes  of  motion.  The  az’muthal  or 
rotational  mode  is  of  inter*;;  here.  In  order  to  predict  the  nature  of 
these  torques  in  quantitative  tirms,  one  must  make  a  thorough  survey  of 
the  highly  complex  wake. 

In  this  report,  we  present  a  qualitative  description  of  the 
wake  behind  a  bluff  body  (like  a  sphere).  Our  remarks  are  drawn  from 
our  own  experiments  as  well  as  the  experiments  of  others. 

The  torque  may  be  estimated  if  the  motion  and  the  dynamic 
parameters  (inertia,  damping,  spring  constant)  of  the  system  are  known. 

In  addition  to  the  fabric  and  hardware  inertia,  the  boundary  layers  con¬ 
tribute  appreciably  to  the  ir.srtia  of  the  system.  The  gases  on  either 
side  of  the  balloon  exert  retarding  torques  on  the  balloon,  the  effect 
being  that  of  a  damping  aerodynamic  torque  and  an  apparent  additional 
moment  of  inertia.  Experimental  values  of  the  damping  aerodynamic 
torque  on  a  sphere  rotating  with  constant  velocity  were  presented  in 
Technical  Report  No.  1  (Ref..  1).  Since  then,  we  have  made  finer  and 
more  complete  measurements ,  and  the  final  results  are  presented  in  this 
report.;  Also,  the  experimental  verification  of  the  formula  of  the  ap¬ 
parent  additional  moment  of  inertia  is  presented  here. 

With  these  aerodynamic  parameters  specified,  a  dynamic  analysis 
with  numerical  examples  is  made  of  the  coupled  rotations  of  the  full 
scale  balloon-gondola  system..  The  magnitude  of  the  forr'r.g  totMue  is 
estimated  for  known  rotations. 
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III.  THE  WAKE  OF  A  MOVING  BhUFF  BODY 


The  flow  of  «  viscous  fluid  around  5  bluff  body  can  be  de¬ 
scribed  briefly  in  the  following  manner.  Viscosity  demands  that  the 
fluid  touching  the  surface  of  the  body  be  stationary  vith  respect  to  the 
body  (no-alip  boundary  condition).  A  small  distance  oway  from  the  sur¬ 
face,  the  flow  is  essentially  potential  (inviccid).  In  other  words,  a 
chin  boundary  layar  covers  the  body,  through  wnich  the  velocity  of  the 
fluid  changes  v^ry  rapidly  from  zero  at  the  surface  to  the  value  of  the 
potential  flow.  The  thickness  of  the  boundary  layer  decreases  with  in- 
cresslng  Reynolds  number  (Re)  Hence,  for  large  values  of  Re,  the 
boundary  layar  is  thin. 

The  strong  viscous  fo'ces  in  end  el  vug  the  boundary  layer  are 
balanced  by  che  pressure  grtlletc  in  the  same  direction  of  the  potential 
flow  outaide  the  layer.  AC  some  point  behind  the  front  of  the  body, 
this  pressure  gradient  may  not  bs  cble  to  cancel  the  viscous  forces  of 
the  boundary  layer.  There,  the  boundary  layer  detaches  from  the  body. 
There  la,  now,  a  surface  starting;  from  the  body  and  extending  behind  it 
in  the  fluid,  on  which  the  velocity  of  Che  fluid  la  equal  to  zaro.  In 
the  two  regions  separated  by  this  surface,  th  flow  is  in  opposite  di¬ 
rections;  i.e.,  we  have  separation  of  che  flo-.  The  velocity  gradients 
normal  to  this  surface  ere  very  high.  Therefore,  this  surface  la  like 
a  vortex  sheet  and,  being  Inherently  unstable  curls  upon  Itself  cre¬ 
ating  a  region  of  highly  concentrated  vorticity;  i.e.,  e  bound  vortex. 
The  vortex  is  fed  from  the  vorticity  generated  at  the  surface  of  the 
body  and,  if  tha  rate  of  generation  of  vorticity  is  greater  than  che 
race  of  its  diffusion  through  the  fluid,  the  vortex  grows  in  strength 
continuously.  Ultimately,  tils  bound  vortex  escapee  and  travels  in  the 
wake,  and  che  cycle  is  repee:ad.  The  overall  result  is  an  occasional 
(often  periodic)  shedding  of  vortices  in  the  vake.  In  svaoary,  the 
shedding  of  vortices  in  the  /eke  Is  due  to  the  >ep*r,.tion  of  the  flow 
and  the  lnsteblllty  of  the  resulting  vortex  sheet  aud  bound  vorticea. 

These  phenomena  have  been  observed  experimentally.  Por  a  cyl¬ 
inder,  the  shed  vortices  are  line  vortices  with  their  axis  parallel  to 
the  axle  of  the  cylinder.  They  ere  ehed  in  tvo  trails.  Since  vortex 
lines  must  be  closed  (like  electric  current  lines),  the  vorticity  of  one 
vortex  in  one  trail  Is  equal  and  opposite  to  ihf  vorticity  of  one  vortex 
in  the  other  trail.  The  two  vortices  meet  at  tie  region  of  the  wake 
corresponding  to  the  edges  o?  the  cylinder  ,  th  us,  forming  e  '•lored  vor- 
tlcity  circuit.  In  the  idea*  cast,  the  tvo  trufls  ere  parallel  and 
neriodlc  constituting  the  so  celled  "Kerman  vertex  street".  These  vor¬ 
tices  exert  forces  on  the  cy  inder,  which  can  cause  lateral  motions  of 
the  cylinder  as  well  as  a  rc  atlon  about  its  exit.  Since  the  vort.clty 
in  one  trail  is  equal  and  op  oeite  to  the  vorticity  in  the  other  trail, 
lunmm  fcrcss  ere  alternating  end  the  resulting  'motions  oscillatory,  rha 
frequency  of  the  oscillation  la  equal  to  tha  l rcquency  of  tha  shedding 
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of  vortices  in  the  wake.  A  curve  has  been  obtained  experimentally,  and 
relates  the  dependence  of  the  Strouhal  number,  St  (a  dimensionless 
parameter  proportional  to  the  frequency),  on  Re  (see  Fig.,  149  of  Ref.  2). 
Indeed,  ic  has  been  shown  that  the  frequency,  in  which  long  wires  "sing" 
under  the  influence  of  the  wind,  can  be  predicted  from  this  curve. 

This  curve  shows  that  St  increases  with  Re  for  100  <  Re  <  700, 
St  is  essentially  constant  for  700  <Re  <  105,  and  St  increases  very 
rapidly  with  Re  for  10^  <  Re  <  10^.  Hence,  the  wake  has  a  rather  well- 
defined  periodicity  in  the  middle  range  of  Re.  In  the  last  range  of  Re, 
small  distrubances  of  the  ambient  flow  will  change  the  wake  radically; 
the  periodicity  of  the  wak-  lost,  and  the  wake  is  very  confused  and 
turbulent . 

The  wake  of  a  bluff  body  with  a  spherical  shape  will  have  the 
same  basic  features  as  that  of  a  cylinder  with  some  differences  in  detail 
due  to  tiu-  differences  in  geometry.  Unfortunately,  the  wake  of  a  sphere 
has  not  been  studied  to  the  same  extent  as  the  wake  of  a  cylinder.  Some 
revealing  experimental  studies  have  been  published  only  recently.  At 
the  present  time,  it  appears  that  there  is  an  adequate  knowledge  of  the 
wake  of  a  sphere  for  Re  up  to  about  1,000  only. 

Magarvey,  _et  al.,  (Refs.  3,  4  and  5)  have  taken  some  excellent 
photographs  of  the  wake  of  a  colored  liquid  drop  ft  ling  freely  in  water. 
Initially,  the  drop  was  accelerating,  but  it  reached  a  constant  velocity 
ultimately.  Figure  l  shows  the  formation  of  the  wake.  As  the  drop  was 
accelerating,  the  flow  began  to  separate  and  a  cylindrical  vortex  sheet 
appeared  (a).  This  sheet,  being  inherently  unstable,  curled  upon  itself 
more  and  more  (b  and  c)  forming  a  bound  vortex  ring.  As  the  vortlclty 
of  the  vortex  ring  Increased  ,  It  began  to  break  off  from  the  drop.  This 
is  shown  in  photograph  (d)  taken  0.01  second  before  the  first  vortex 
ring  was  discharged  in  the  wake.  Note  that  the  ring  is  discharged  at  a 
distance  behind  the  drop  of  the  order  of  magnitude  of  the  diameter  of 
the  drop.  F<gure  2  shows  the  wake  of  a  drop  falling  with  constant  ve¬ 
locity  corresponding  to  Re  =  360  ((b)  is  10  seconds  later  than  (a)  and 
shews  the  diffusion  of  the  vortex  elements)..  The  vortex  elements  are 
distorted  rings  interconnected  with  line  vortices.  There  is  a  definite 
periodicity  and  a  separation  of  the  vortices  in  two  trails.  The  vortlc- 
ity  of  the  rings  in  both  trails  is  in  the  same  direction.: 

Magarvey  obtained  an  experimental  curve  for  the  dependence  of 
St  on  Re  for  270  <  Re  <  390.  As  Re  was  increased  further,  the  wake  ut- 
gan  to  lose  its  periodicity.  This  is  shown  in  Figure  3  io<  »<•  1040 

Figure  4  shows  representative  results  froai  our  own  rough  ex¬ 
periments.  A  ping-pong  ball  was  placed  in  a  uniform  stream,  vh  speed 
was  varied  so  that  values  of  Re  up  to  6,000  were  obtained.  Ink,  flowing 
out  from  small  h'les  on  the  equator  of  the  ball,  smde  the  wake  ..ole. 

At  Re  >  600  (Fig.  4a),  the  wake  was  similar  to  chat  obtained  by  Magarvey, 


FIGURE  2 


THE  WAKE  OF  A  LIQUID  DROP  OF  REYNOLDS  NUMBER  360. 
(COURTESY  OF  "THE  PHYSICS  OF  FLUIDS") 


FIGURE  3 


THE  WAKE  OF  A  LIQUID  DROP  OF  REYNOLDS  NUMBER  HW! 
(COURTESY  OF  "CANADIAN  JOURNAL  OF  PHYSICS"! 


and  it  had  a  definite  periodicity.  At  Re  =  1200  IFig.  4b),  some  peri¬ 
odicity  was  still  visible.  At  Re  =  2400  (Fig.  4c),  the  wake  was  quite 
confused  with  some  regions  of  concentrated  vorticity.  At  Re  =  4800 
(Fig.  4d) ,  the  wake  was  almost  uniformly  turbulent  and  all  periodicity 
disappeared. 

This  apparent  disorder  in  the  wake  of  a  sphere  at  values  of 
Re  above  1,000  must  be  due  to  the  presence  of  a  high  degree  ol  ambient 
turbulence,  or  a  very  sensitive  dependence  of  St  on  Re,  or  both.  Accu¬ 
rate  and  well-controlled  experiments  will  clear  up  this  point.;  At  the 
present  time  one  thing  is  rather  clear;  i.e.,  for  Re  in  tne  range  of 
about  one  to  a  few  thousands,  the  wake  of  the  ephere  is  incipiently  un¬ 
stable,  and  very  small  distui nances  of  the  uniform  stream  (like  Inherent 
turbulence,  small  ambient  flow  disturbances,  etc.)  can  alter  the  form 
of  the  wake  and  destroy  its  periodicity.  Inis  has  been  confirmed  by 
other  investigators  (verbsl  communications).; 

As  in  the  case  of  the  cylinder,  the  wake  of  the  sphere  exerts 
forces  and  torques  on  the  sphere.  Therefore,  if  the  sphere  is  free,  it 
can  assume  quite  a  few  modes  of  lateral  motion  and  rotation.  The  same 
phenomena  will  occur  for  a  body  liks  a  balloon.  These  facts  have  been 
confirmed  by  our  experiments  with  balls  and  scaled  balloons  falling 
freely  In  a  water  tank  (Ref.  1). 
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IV.  DAMPING  AERODYNAMIC  TORQUE  FOR  A  ROTATING  SPHERE 


Since  Che  azimuthal  rotation  of  balloons  is  especially  crucial 
while  they  are  at  or  near  their  floating  altitude,  we  have  restricted 
our  attention  to  the  case  of  fully  inflated  balloons.  A  balloon  at 
floating  altitude  is  dealized  as  a  spherical  shell  with  an  extremely 
boall  thickness  to  radius  ratio.  The  spherical  shape  is  believed  to  be 
close  enough  to  the  cnion-shape  of  actual  balloons  to  make  the  results 
applicable,  while  preserving  geometric  simplicity  for  ease  of  analysis. 
It  is  further  assumed  that  damping  is  independent  of  the  small  values  of 
vertical  velocity. 

A  simple  experiment  (Ref.  1,  p.  36)  was  previously  described 
for  obtaining  tbe  torque  required  to  overcome  viscous  skin  friction  on 
a  steadily  rotating  sphere.  The  goal  of  this  experiment  was  to  provide 
a  basis  for  estimating  the  damping  aerodynamic  torque  which  opposes  the 
rotations  of  rising  and  falling  balloons. 

The  results  of  the  first  experiment  were  sufficiently  encour¬ 
aging  to  wartant  refining  the  experiment.  While  no  better  means  of 
measuring  small  steady  torques  was  found,  more  precise  measurements  have 
been  carried  out  on  a  greater  number  of  models  which  were  more  highly 
spherical.  The  data  extend  over  a  wider  range  of  variables  and  are  more 
reliable  than  the  previous  data. 


A.  EXPERIMENTAL  SET-UP 


A  stainleas  steel  cylindrical  tank--3  feet  in  diameter, 

2.5  feet  deep--was  filled  with  water.  Brookfield  viscometers  (Model  RVF) 
were  used  to  spin  the  models  at  constant  speed  about  a  vertical  axis  and 
to  indicate  the  torque  required  to  maintain  the  steady  rotation.;  One 
viscometer  had  operating  speeds  of  ’00,  50,  20,  and  10  rpm;  the  other, 

20,  10,  4,  and  2  rpm.  Each  in  turn  was  carefully  mounted  with  its  shaft 
at  the  center  of  the  tank.  The  models  were  completely  filled  with  water 
and  submerged  so  that  tho  "north  pole"  was  about  3  inches  beneath  toe 
free  surface. 

Eccentricity  was  minimized  by  very  careful  mounting.-  still  a 
small  AC  disturbance,  due  to  eccentricity,  appeared  in  some  viscometer 
readings.  The  average  readings  over  a  cycle  should  not  he  effected  by 
it. 


Sufficient  time  was  allowed  after  starting  tht  synchronous 
motor  in  the  viscometer  for  the  starting  transients  to  damp  o-» .  (see 
further  discussion  below) . 
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B.  SPHERICAL  MODELS 


The  models  were  spherical  aluminum  floats  (Fig.,  5)  .  Each 
float  consisted  of  two  hemispheres,  of  about  0.025  inch  wall  thickness, 
welded  together  at  the  equator.  The  welded  seams  were  filed  down  flush, 
and  each  sphere  was  sprayed  with  clear  lacquer. 

The  nominal  diameters  were  2,  3,  4,  5,  6,  7,  and  B  inches,  but 
the  two  largest  spheres  could  not  be  used.  They  put  enough  weight  on 
the  viscometer  bearing  to  increase  bearing  drag  prohibitively.  The  ac¬ 
tual  equatorial  diameters  as  measured  and  the  mass  parameters  (for 
future  reference)  are  tabulated  in  Table  1.  Generally,  the  equatorial 
and  polar  diameters  agreed  w'  in  about  1/64  inch,  but  the  polar  diameter 
of  Sphere  Number  5  was  about  3/64  inch  leas  than  its  average  equatorial 
diameter. 


TABLE  I 

DIAMETERS  AND  MASS  PARAMETERS  OF  SrHERICAL  FLOATS* 

Moment  of  Inertia 


Sphere  No. 

Diameter  (CM) 

Maas  (GM) 

(GMCM2)** 

1 

5.02 

19.5 

81.9 

2 

7.58 

67.3 

646 

3 

10.12 

128.0 

2,190 

4 

12.66 

197.6 

5,290 

5 

15.16 

186.5 

7,170 

6 

17.66 

399.7 

20  800 

7 

20.20 

478.1 

32,500 

‘Obtained 

from  the  Chicago  Float 

Works  2330 

South  Western  Avenue, 

Chicago  8,  Illinois. 

“Computed  from  Moment  of  Inertia  *  1/6  (Mass) (Diameter)2 . 
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EXPECTED  FORM  OF  THE  DAMPING  TORQUE 


In  Technical  Report  No.  1  (Ref.,  1)  the  form  of  the  aerodynamic 
damping  torque  on  a  rotating  sphere  was  deduced  from  theoretical  consid¬ 
erations  of  the  boundary  layer  uea.*  a  flat  plate  moving  with  constant 
velocity.  Here,  the  form  r2  the  torque  will  be  derived  from  the  flow  of 
a  rotating  dish,  which  has  many  similarities  with  the  flow  of  r  rotating 
sphere  and  has  been  completely  solved.. 

Figure  6  shows  the  streamlines  of  the  flow  in  a  meridional 
plane  of  an  infinite  disk  and  a  sphere  of  radius  a,  both  rotating  with  a 
constant  angular  velocity  to.  except  in  the  neighborhood  of  the  equator, 
the  flow  of  the  sphere  is  quite  similar  to  the  flow  of  the  disk.  Through¬ 
out  the  disk  the  boundary  layer  thickness  is  constant.  For  the  sphere, 
the  two  boundary  layers,  which  start,  from  the  poles,  grow  towards  the 
equator,  where  they  collide  and  erupt.:  (Nor  is  the  flow  of  the  rotating 
sphere  near  the  equator  like  the  flow  of  a  rotating  cylinder,  since  in 
the  latter  there  is  no  secondary  outflow  at  all.) 

For  the  disk,  the  angular  velocity  of  the  fluid  decreases  to 
about  0.03u  in  a  distance  $  away  from  the  disk,  which  is  given  by  (see 
Ref.;  2,  Section  43) : 


s  = 


(1) 


where  ft  and  p  are  the  viscosity  and  density  of  the  fluid,  respectively. 
Taking  S  as  the  boundary  layer  thickness  for  the  sphere  and  using  the 
local  velocity  v  =  wasinO,  we  find  that  the  shearing  stress  (*  U  v/  4  ) 
acting  on  the  aphere  at  0  ia  approximately  equal  to  0.222  (ppV'  2toi 2 
aslnO.  Then  the  damping  torque,  D,  which  is  equal  to  the  summation  over 
the  sphere  of  the  product  of  the  shearing  stress  by  the  area  by  the  dis¬ 
tance  fro.-,  *.h.  axis,  ia  approximately  given  by:' 

Jrir 

Td. 222  ( pft)l/ 2  u> 3/2  »  sin  ?]  (2  V  a2  sin  8  d8)<»  sin  8). 
o  L  J 

Therefore, 

D  *  1.9  (pp)VZ  CO3'2  a'+.  (2) 


Let  us  define  the  pertinent  Reynold's  number,  R,  by; 


R 


i3) 
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FIGURE  6 


STREAMLINES  IN  A  MERIDIONAL  PLANE  OF  (a)  ROTATING  DISK  AND  (b) 
ROTATING  SPHERE 


I 

J 

! 


and  a  dimensionless  drag  coefficient:,  CQ,  by:; 


D 


P  *»2  cn 


*4) 


Then  Equation  2  can  be  put  in  the  form:' 


1.9  R 


-1/2 


(5) 


The  form  of  Equations  2  and  S  will  be  the  guide  in  correlating 
the  experimental  results . 


0.  EXPERIMENTAL  RESULTS  AND  DISCUSSION 


The  experimental  data  are  tabulated  in  Table  II.  Each  datum 
point  tabulated  la  the  average  of  at  least  three  tests.  The  repeatibil- 
ity  of  the  torque  measurements  was  about  1  3.6  dyne-cm. 

Figure  7  shows  the  experimental  data  plotted  non-dimens  Iona lly 
as  C0  versus  R.  In  the  experimenta  R  was  varied  from  670  to  27,300;  the 
torque  varied  from  10  dyne-cm  to  7130  dyr.e-cm  (i.e.,  by  an  amount  ap¬ 
proaching  three  orders  of  magnitude).  Over  this  wide  range,  the  data 
are  very  well  described  by  the  following  simple  law: 


=  3.6  R 


-1/2 


(6) 


which  appears  as  the  solid  line  in  Figure  .  Combined  with  Equation  4,. 
this  leads  to  the  formula  for  the  data:- 


D 


3.6 


1/2 


4  3/2 

a  w 


(7) 


This  compares  well  with  the  scanty  previous  data  suBnarized  by  Equation 
4.3  of  Reference  1  (p.  38).  As  the  spheres  were  filled  wltn  water, 
Equation  7  represents  the  total  torque  frtw  externa!  internal  flows. 

Equation  6  shows  that  C0  is  a  function  of  the  Reynold's  rusher 
ouly.  This  is  expected  since  the  drag  here  is  entirely  viscous.  TK- 
nonlinear  dependence  of  the  torque  on  u  (D  proportional  to  instead 
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TABLE  II 


SunivuvV  Or  EXFa&InfcWIAL  UA1A 


Sphere 

No. 

Ui 

rad/sec 

Torque,D, 

dvne-cm 

D 

C _ 0__ 

,A  .1/2  3/2  4 

kPu)  Ui  a 

CD  -  — 1"5 

PiO  a 

1 

10.48 

471 

3.54 

0.0431 

,  .....„3 
o  .  /  i  iu 

5.24 

162 

3.44 

0.0593 

3.36-10? 

2.096 

37.7 

3.16 

0.0864 

1.34-10-’ 

1.048 

10.06 

2.38 

0 .0919 

6  7l'l0J 

2 

10.48 

2440 

3.52 

0.0284 

1.54' 104 

5.24 

87? 

3.58 

0.0409 

7  .68  •  103 

2.096 

220 

3.55 

0.0640 

3.08-103 

1.048 

73.6 

3.36 

0.0857 

1.54-103 

3 

10.43 

7130 

3.23 

0.0196 

2.73-104 

5.24 

2945 

3.78 

0.0324 

1 .37 • 104 

2.096 

733 

3.72 

0.0504 

5.46-103 

1.048 

253.5 

3.63 

0.0695 

2.73' 101 

4 

2.096 

1746 

3.62 

0.0392 

8.55-103 

1.048 

675 

3.96 

0.0606 

4.27-103 

5.24 

6150 

3.43 

0.0222 

2.40-104* 

2.096 

1680 

3.70 

0,03/9 

9 .56 • 103* 

1.048 

589 

3.67 

0.0531 

4.78-103* 

5 

2.096 

3650 

3.91 

0.0334 

1 .37-104* 

1 .048 

1210 

3.67 

0.0443 

6.26-103* 

0.4192 

331 

3.97 

0.0758 

2.74103* 

2.096 

4020 

4  07 

0.0367 

1.23-1C4 

1.048 

1393 

3.98 

0.0508 

6.13- 103 

0.4192 

392 

4.43 

0.0895 

2.45-103 

*  For  entries  with  end  without  asterisk,  the  water  temperature  was  26°C 
and  21°C,  respectively.  The  resulting  change  in  viscosity  caj  id 
slightly  different  values  of  R  for  the  same  values  of  a  and  Ui- 
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AMPINO  TORQTF  COEFFICIENT  (C  )  VERSl  S  ROTATIONAL,  REYNOLDS  NCMHER  (R) 


of  ui  as,  fui  instance,  in  the  case  of  a  rotating  cylinder)  is,  of 
course,  due  to  the  existence  of  a  boundary  layer  whose  thickness  is 
proportional  to  to'1'  (see  Section  IV-C) . 

Thus,  it  has  been  verified  that  the  dependence  of  U  on  the 
various  parameters  of  the  problem  is  of  the  form  anticipated  from  theo¬ 
retical  considerations  in  Section  IV-C.  The  constant  of  proportionality 
has  been  determined,  and  it  is  essentially  independent  of  the  Reynold's 
number  for  the  wide  range  of  Reynold's  numoers  considered  here  (from 
670  to  27,300).  Considering  the  fact  that  both  the  internal  and  the 
external  flows  have  contributed  equally  to  the  torque  given  oy  Equation 
7  (see  Section  IV-E),  it  ’  seen  that  there  is  a  satisfactory  agreement 
between  theory  and  experiment  even  in  the  numerical  value  of  the  con¬ 
stant  of  proportionality  (1.9  by  theory  and  3.6/2  =  i.8  by  experiment). 
But  this  should  rather  be  considered  as  a  coincidence,  since  in  tne 
theory  the  value  of  this  constant  depends  on  what  exactly  is  taken  as 
the  thickness  of  the  boundary  layer.  In  Section  IV-C  we  defined  the 
thickness  of  the  boundary  layer  as  being  equal  to  the  distance  from  the 
disk  at  which  the  angular  velocity  of  the  fluid  decreases  to  0.03  u>  . 

We  can  define  this  thickness  at  a  different  velocity  decrease  (say 
0.1  a)  or  0.0’  a>) ,  and  the  value  of  the  constant  in  Equation  2  will  be 
different  but  the  form  of  this  equation  will  be  the  same.  All  that  the 
present  theory  is  expected  to  give  is  the  correct  order  of  magnitude 
for  the  constant  of  proportionality. 


E.  INTERNAL  FLOW 

As  mentioned  at  the  beginning  of  Section  IV-A,  the  models  were 
completely  filled  with  water.  Therefore,  there  was  an  internal  as  well 
as  an  external  flow.  The  two  flows  can  differ  greatly  in  the  following 
respect.  In  the  outer  flow,  the  effects  of  viscosity  continuously  dif¬ 
fuse  outward  so  that  the  fluid  is  always  sheared  at  the  surface.  How¬ 
ever,  Inside  the  sphere,  the  effects  of  viscosity  accumulate,  and  if  the 
steady  sphere  rotation  persists  long  enough,  then  the  fluid  inside  may 
rotate  with  the  sphere  as  a  rigid  body.  When  this  happens,  there  is  no 
shearing  at  the  inner  surface  of  the  sphere  and,  therefore,  no  torque 
acting  on  it.; 

The  effects  of  viscosity  diffuse  inward  to  a  depth  equal  to 
the  radius  a  in  a  time  f  of  the  order  pallet.  It  will  be  at  least  of 
this  order  of  time  before  rigid  body  rotation  of  the  internal  fluid  is 
possible.  For  the  smallest  sphere  used  in  the  cr. perimer.es ,  X  -s  of 
the  order  of  11  minutes.  All  readings  were  taken  in  much  snorter  times 
aits'  starting  the  rotation,  so  we  infer  that  the  internal  boundary 
layer  was  relatively  thin  and  the  data  contain  contributions  '  in¬ 
ternal  end  external  zlows. 
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In  order  Co  apply  the  pretent  result*  to  balloons,  which  have 
different  fluids  Inside  and  outside,  it  it  necessary  to  determine  the 
relative  contributions  of  internal  and  external  flows  to  the  total 
torque.  An  attempt  to  perform  this  separation  of  effects  experimentally 
was  not  successful.  The  models  were  first  filled  with  water  and  rotated 
in  air;  then  filled  with  air  and  rotated  in  water.  In  both  cases,  there 
was  Coo  such  verticil  lead  or.  the  viscometer  bearing  to  yield  any  re¬ 
sults  of  value.  This  difficulty  could  probably  be  avoided  in  the  future 
by  using  alcohol  and  water  for  the  two  working  fluids,  and  this  experi¬ 
ment  is  contemplated. 

On  the  other  hand,  'r.cs  in  the  experiment  both  internal  and 
external  boundary  layers  were  very  thin,  it  seems  reasonable  to  assume 
that  the  internal  and  external  flows  contributed  the  same  aaount  to  the 
total  measured  torque.  This  assumption  is  supported  by  the  results  of 
Lamb  (kef.  6),  who  analyzed  theoretically  the  ssiall  amplitude  rotational 
oscillations  of  spherical  shells  with  fluid  either  inside  or  outside  the 
shell.  He  found  that  the  hydrodynamic  torques  from  internal  and  external 
flows  are  identical  when  the  boundary  layers  are  very  thin. 
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APPARENT  ADD IT  TONAL  MOMENT  OF  I!  'ERTIA  FOR  A  ROTATING  SPHERE 


The  preceding  section  has  de*' t  with  the  torque  arising  from 
the  steady  rotation  of  the  zero  position.  However,  an  additional  reac¬ 
tive  torque  is  associated  with  rotational  accelerations,  which  Is  the 
subject  of  this  section. 

When  any  solid  body  moves  with  constant  velocity  in  an  incom¬ 
pressible  and  inviscid  fluid,  the  net  force  exerted  by  the  fluid  on  the 
body  is  equal  to  zero.  On  the  other  hard,  when  the  body  stoves  with  a 
constant  acceleration,  the  surrounding  fluid  exerts  on  the  body  a  net 
retarding  force  which  is  proportional  tc  "he  acceleration.  This  effect 
has  become  known  as  th  *  "apparent  or  induced  mass"  effect,  since  the 
body  appears  to  have  more  mss  when  aec&lerated  in  the  fluid. 

For  a  smooth  sphere  accelarat lag  angularly  about  its  axis  in 
a  viscous  fluid,  a  similar  type  of  affect  occurs.  However,  the  physical 
mechanism  causing  this  affact  is  now  different.  Here,  viscosity  demands 
a  "no-slip"  boundary  c  ndition  on  the  surface  of  the  sphere  (or  balloon). 
Consequently,  whenever  the  balloon  Accelerates  angularly,  some  of  the 
surrounding  air  and  contained  helium  mus  r.  accelerate  with  the  balloon. 
Hence,  a  net  retarding  torque  is  exerted  on  the  balloon  by  the  air  and 
helium.  In  addition  to  a  damping  effect,  this  torque  contains  an  in¬ 
ertial  effect  referred  to  as  the  "apparent  additional  moment  of  inertia" 
of  the  balloon. 

For  most  solid  bodies,  with  thin  boundary  layers,  this  increase 
in  Inertia  is  negllgib’e  compared  to  that  already  in  the  system  by  vir¬ 
tue  of  solid  mss.  However,  fully  inf  la  fed  balloons  have  a  very  sull 
moment  of  inertia  (for  bodies  of  their  size),  and  the  apparent  additional 
moment  of  inertia  is  sgnlf leant.  The  fallowing  experiment  was  intended 
tc  provide  an  estimate  of  the  effect. 


A,  DESCRIPTION  OF  THE  "XPERIWNT 

In  the  experlMnt,  a  sphere,  i amerced  in  water,  was  given  a 
constant  initial  rotation  and  than  tha  d  scay  of  tha  rotation,  due  to  the 
viscosity  of  water,  wa;  recorded.  From  this  record,  the  apparent  addi¬ 
tional  moment  of  inert? a  of  tha  sphara  vis  calculated. 

A  photograph  of  tha  axparlawntal  set-up  Is  chown  in  Figure  6. 

A  tank,  3  feet  in  diameter,  was  filled  with  clean  tap  water  to  a  level 
of  2.5  feat.  Soma  of  :he  spherical  mode’s,  described  in  Section  IV-B, 
were  threaded  on  tha  end  of  a  12-inch  lea?  and  1/4-lnch  dlamate  pre¬ 
cision  ground  stainless  steal  shaft.  Th:  jhaft  was  placed  in  a  4 -inch 
long  stainless  steal  cylindrical  bearing  <■  older  having  two  precis. ./a 
stainless  steel  ocarinas,  ons  at  aach  an:  This  system  was  suepandad 


3rttur  B.litrtr.  Jar. 


Arthur  D.'D.iUlr  Jhn 


21 


vertically  along  the  centerline  of  the  water  tank  by  mounting  the  bear* 
ing  holder  to  a  heavy  wood  framework  (see  photograph1' .  The  upper  bear¬ 
ing  was  a  thruat  bearing  supporting  the  net  weight  of  the  model  and 
•haft,  while  the  lower  was  a  simple  radial  bearing  to  minimize  lateral 
motions.  From  the  manufacturer's  data,  the  frictional  torque  of  the 
bearing  is  about  200  cm  dynes  for  a  thrust  of  200  gm.  The  level  of  the 
water  in  the  tank  was  about  5  Inches  above  the  top  of  the  model.. 

The  upper  end  of  the  shaft  protruded  upwards  from  the  bearing 
holder  about  2  Inches.  There,  a  thin  aluminum  circular  disk  (to  be  de¬ 
scribed  later)  was  attached  to  the  shaft  In  a  horizontal  position  and, 
above  the  disk,  a  tee-shaded  driving  yoke.  Located  above  the  yoke  was 
a  driving  mechanism  consisting  of  a  laboratory  stirrer  with  variable 
speed  coupled  to  a  "Zero-Max"  speed  reducer.  The  output  of  the  driving 
mechanism  was  variable  from  zero  to  about  three  revolutions  per  second.; 
On  the  output  of  the  speed  reducer,  there  was  a  simple  two-pin  clutch 
which  engaged  the  tee-shaped  driving  yoke.  After  tne  system  was  brought 
up  to  a  preselected  speed,  the  clutch  was  disengaged  and  the  rotating 
mass  was  free  to  slow  down.  The  decay  of  tne  rotation  of  the  spherical 
model,  due  mostly  to  viscous  effects,  was  recorded  by  an  optical  system. 

The  light  source  (prefocused  beam)  and  the  detector  (photo- 
duo-diode)  of  the  optical  system  are  mounted  opposite  to  each  other  in  a 
C-shaped  frame.  In  the  experiment,  this  frame  was  placed  so  that  the 
aluminum  disk  interrupted  the  light  path.  The  disk  had  small  holes, 
placed  evenly  on  a  circumference,  through  which  the  light  beam  passed  as 
the  disk  rotated.  Thus,  a  pulsed  signal  was  generated  by  the  detector, 
which  was  amplified  and  displayed  by  a  time-base  event  recorder.  The  re¬ 
sulting  data,  consisting  of  a  sequence  of  sharp  pulses,  gave  the  angle 
of  rotation  versus  time.  From  these,  the  angular  velocity  of  rotation 
was  calculated  as  a  function  of  time. 

Wobble,  due  to  eccentricity,  was  minimized  oy  very  careful 
mounting  and  by  inserting  soft  rubber  in  the  mounting  of  the  bearing 
holder  to  the  wood  framework.  To  reduce  bearing  loads,  the  models  were 
filled  with  vater.  All  experiments  were  performed  in  calm  water. 

The  initial  speeds  *era  selected  so  that,  foi  ihe  most  part  of 
the  duration  of  the  decay,  the  Reynold's  number  was  within  the  range  xn 
which  the  damping  aerodynamic  torque  was  measured  in  Section  IV  (700  to 
27,000) . 


B.  gXFlCTgg  FORM  OF  TH8  APPARENT  ADDITIOHAL  MOigMT  OF  1HKRT1A  AND 
OTHER  THEORETICAL  COWS1DERATIOMS 


Consider  the  sphere  rotating  with  a  constant  angulax  'elocity 
lit.  It  has  been  pointed  out  in  Section  IV-C  that  tne  flow  1«  <’  ce 
•  teller  to  that  of  a  rotating  disk  (.see  Fig.  6).  Taking  the  tni:knees 
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of  the  boundary  layer  &  as  given  by  Equation  1,  the  total  angular  mo¬ 
mentum  of  the  fluid  in  the  boundary  layer  of  the  sphere  is  given  ap¬ 
proximately  by: 

I  (/>)(“^)(*2  sin2  S)  ( 2  rr  c}  sin  6  d  0)  =19  L}1' C  e!* 

Jo 


The  torque  exerted  by  the  fluid  on  the  sphere  is  the  steady  torque  l) 
dealt  with  in  the  preceding  section.; 

Suppose,  now,  that  the  angular  velocity  of  the  sphere  is 
changed  by  an  amount  ^(*>  in  a  time  fat.  The  surrounding  viscous  fluid 
will  resist  this  change,  and  in  addition  to  the  steady  torque  0,  it  will 
exert  on  the  sphere  a  torque  equal  to  the  time  rate  of  change  of  its 
angular  momentum  as  found  above.  Therefore,  the  apparent  additional 
moment  of  Inertia  iA  of  the  sphere  is  given  aDp'-oximately  by: 


«  19  (,A)l/2a>-l/2.4 


:*> 


Actually,  this  is  a  quasi-steady  state  analysis,  not  always 
valid.  It  has  been  assumed  that  the  time  fat,  in  which  the  change  fat,} 
occurs,  is  much  larger  than  the  time  required  for  a  disturbance  occur¬ 
ring  at  the  surface  of  the  sphere  to  diffuse  through  the  boundary  layer, 
so  that  the  boundary  layer  remairs  just  about  fully  grown  during  the 
change  in  to  ,  This  latter  time  is  of  order  p  /  fl  -  1/a )-  The  condi¬ 
tion  for  the  validity  of  Lhe  abo"e  analysis  is  a  condition  on  the  magni¬ 
tude  of  the  angular  acceleration  of  the  sphere,  which  can  be  stated  as; 


A  a) 
A  i 


<  <  U) 


(9) 


Neglecting  bearing  friction,  the  rotation  of  the  sphere  must 
satisfy  tne  following  differential  equation; 


_d 

dt 


r~ 

[<I  *  2Ia)u> 


2D 


0 
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where  I  is  the  material  moment  of  inertia  of  the  spherical  shell.  The 
factor  of  2  appearing  in  front  of  IA  and  D  accounts  for  both  externa) 
and  internal  flows,  which,  it  is  assumed,  contribute  equally  (see 
Section  lv-E) .  IA  will  be  taken  in  the  form  of  Equation  8  but  with  an 
unknown  constant  of  proportionality  to  be  determined  from  the  experiment. 
D  will  be  taken  as  found  in  Section  IV.  That  is:: 


,  „  .  .1/2  4  -1/2 

IA  =  K  s/9/2)  a  u) 


(11) 


1.8  a4  U>3/2 


(12) 


Let  u>  denote  the  value  oi  u  at  t  =  0.  Substituting 
Equations  11  and  12  in  Equation  10  and  introducing  the  following  dimen¬ 
sionless  quantities:- 


i  a  „4  „  1/2 

3.6  a  u)Q 


(13) 


,  ,1/2  T 
<^0>  1 


(14) 


the  solution  of  Equation  10  is:- 


&f  ■  '«  ‘  [|‘  *«>2 


(13) 


Notice  that  the  parameter  ot  is  the  ratio  of  cue  material  mo¬ 
ment  of  inertia  of  the  spherical  shell  to  the  apparent  additional  moment 
of  inertia  due  to  the  external  or  internal  flow  at  initial  conditions.; 

X  is  a  dimensionless  time..  The  experimental  results  will  be  plotted  as 
(u0/w  )1/2  versus  t  .,  Then  from  the  family  of  curves  giv:.  Ly  Equation 
13,  the  cuive  that  fits  the  experimental  results  best  wi'l  be  selected.; 
This  determines  the  value  of  cl  from  which  the  unknown  constant  K  can  be 
computed . 
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C.  EXPERIMENTAL  RESULTS  AND  DISCUSSION 

The  results  of  three  experimental  runs  are  tabulated  in  Table 
III.  The  three  models  were  filled  with  water.;  The  results  satisfy  the 
following  three  conditions,  and  therefore,  the  theory  developed  in  the 
preceding  section  is  applicable:; 

2 

1.  The  total  time  tor  each  run  is  much  less  than  p  &  /  pt 
and,  therefore,  the  assumption,  that  the  hydroaynaraie 
effects  of  the  internal  and  external  flows  are  equal, 
seems  reasonable.; 

2 

2.  The  Reynold's  number  (p  i o  &/u)  is  within  the  range 
(700  to  27,000)  in  which  the  damping  aerodynamic  torque 
has  been  determined.- 

3.  The  condition  of  small  angular  accelerations,  as  stated 
by  Expression  9,  is  satisfied. 

The  experimental  points  are  shown  in  the  normalized  form 
{ilo  o/ ^  versus  X  ]  on  Figures  9-11.  Notice  that  there  is  an 
oscillation  due.  to  a  slight  wobble  of  the  models.  Three  theoretical 
curves  (Equation  15)  are  drawn  for  each  experimental  run.  In  each  case, 
the  middle  curve  represents  the  best  value  of  the  parameter  d.-  The 
value  of  Ot  for  each  run,  as  well  as  the  corresponding  value  of  .tie  un¬ 
known  constant  K  (computed  from  Equation  14),  is  shown  in  Table  IV .;  The 
average  value  of  K  for  the  three  runs  is  32.  Therefore,  for  either  ex¬ 
ternal  or  internal  flow:- 


IA  =  22  a*  u/1/2 


Notice  that  the  theoretical  and  experimental  values  of  the 
constant  of  proportionality  in  the  formula  of  IA  (Equations  8  and  lb) 
agree  quite  well  (see,  also,  the  last  paragraph  of  Section  IV-l>)  .; 

We  believe  that  there  is  enough  evidence,  in  the  above  experi¬ 
mental  results,  supporting  our  contention  that  Equation  lb  will  give,;  at 
least,,  a  fair  estimate  of  IA  for  the  internal  and  external  flows  under 
the  above  stated  three  conditions. 
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FIGURE  11  CORRELATION  OF  EXPERIMENTAL  POINT?  (SPHERE  #5) 
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TABLE  IV 


EXPERIMENTAL  DETERMINATION  OF  THE  CONSTANT  K 


■  X 


Sphere  No. 


txperimentai 
Value  of  Ot 


Computed 
V  a  J  u  e  o  t  K 


Average  Value  of  K>  22 
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VI.  THE  AEKODvHAMIC  PARAMETERS  OF  FULLY  INFLATED  BALLOONS 


It  is  assumed  that  a  fully  inflated  balloon  at  high  altitudes 
can  be  idealized  as  a  spherical  shell. 

For  a  balloon,  the  inside  ana  outside  gace^are  different 
(helium  and  air).  The  ratio  of  the  value  of  (pp  )  ’ '  ‘  for  helium  to 
that  for  air  is  equal  to  about  0.40  and  is  independent  of  the  altitude 
of  the  balloon.  Applying  the  results  of  tne  preceding  two  sections,  we 
obtain > 


0  =  2,5  (pp)1'2  a4  W3/2 

(17) 

t  ,,  .1/2  4  -1/2 

IA  -  31  (pp)  a  ut 

(18) 

where  p  and  p  pertain  to  air  and  a  is  the  equivalent  radius  of  the 
balloon.  Equations  17  and  18  are  representations  of  Equations  7  and  16, 
respectively,  for  a  helium  filled  sphere  in  air. 

Equation  17  is  valid  for  a  steady  rotation  which  satisfies  the 
tol lowing  two  conditions: 

2 

1. :  The  Reynold's  number  (2>  <0  •  //<  )  1,6  in  Lhe  r*nB* 

700  to  27,000. 

2. ;  The  time  of  the  duration  of  the  steady  rotation  must 

be  much  less  than  p  / p  .  (For  times  much  larger  than 
p  a^ / p  ,  the  contribution  of  the  internal  flow  is  negli¬ 
gible,.  and  2.5  in  Equation  17  should  be  replaced  bv  1.6.) 

The  above  .uo  expressions  are  also  applicable  to  an  unsteady 
rotation  which  sat.sfies  the  following  additional  condition  (see 
Sect  lot  IV - B)  . 

3.  The  angular  acceleration  must  be  small,  i.e."> 


dw 

dt 


<  <  U> 


(19) 


Consider  a  balloon  with  a  radius  of  about  80  ieet  near  cs: 1  lug 
at  au  altitude  of  60  000  feet.  At  this  altitude,  the  values  of  the 


1 1 
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kinematic  viscosity  (fi/p)  for  the  atmosphere  and  the  hen 
the  balloon  are  about  3.4  x  10"3  and  2.4  x  10"2  ft2/sec,  re 


Hum  Inside 
respectively . 

Therefore,  condition  2  states  that  the  duration  of  a  steady  rotation 
should  be  much  smaller  than  70  hours,  which  is  easily  satisfied  in 
practice.  Condition  1  requires  that  the  angular  velocity  should  not 
exceed  0.9  rau/sec,  which  is  also  satii cied  in  practice.  As  for  condi¬ 
tion  3,  it  it  satisfied  by  the  rather  smooth  oscillations  of  the 
November  1959  Project  Strato-Lab  Flight  but  not  by  the  violent  oscilla¬ 
tions  of  earlier  <1955  and  1935)  balloon  flights  (see  Figs.  2,  3,  tind 
4  of  Ref .  1)  . 


Table  V  shows  the  c.iues  of  0  and  1^  for  two  values  of  u> , 
representing  roughly  the  maximum  and  minimum  angular  velocity  of  the 
above  Strato-Lab  Flight  at  ceiling.  The  radius  of  the  balloon  was  about 
80  feet.  The  weight  of  the  balloon  fabric,  constructed  of  2  mil  poly¬ 
ethylene,  was  1100  lbs.  Therefore,  the  moment  of  inertia  of  the  material 
of  the  balloon  was  about  1.5  x  10^  slugs  ft2. 


TABLE  V 

VALUES  OF  THE  AERODYNAMIC  EFFECTS  FOR  A  BALLOON  OF 
80  FEET  IN  RADIUS  AT  AN  ALTITUDE  OF  80.000  FEET 


Angular  Velocity 
U)  (rad /min) 

0.2 

0.04 


Damping  Aerodynamic 
Torque-D  (ft  lbs) 

0.09 

0.008 


Apparent  Additional 
Moment  of  Inertia 
IA  (slug,  ft2) 


10 


2.3  x  10J 


Thus,  the  apparent  additional  moment  of  inertia  is  comparable 
to  that  of  the  feurli.  In  the  range  of  angular  velocities  indicated  in 
Table  V.  Considering  the  large  size  of  the  baHoon,  the  total  moment 
of  inertia  and  the  damping  torque  (aerodynamic)  arc.  indeed,  small. 

Small  moment  of  inertia  makes  it  easy  to  set  the  balloon  in  rotation, 
and  small  damping  results  in  a  very  slow  decay  of  the  rotation.  Hence, 
balloons  can  be  very  easily  rotated. 

finally,  we  like  to  present  a  v*iy  brief  review  of  the  exist¬ 
ing  literature  on  tne  aerodynamic  parameters  of  a  rotating  sphe  s  and 
discuss  their  applicability  to  balloon  rotations.  Lamb  (Ref.  Cl  ban 
dealt  with  the  p  oblsm  of  a  sphere  with  sinusoidal  rotation  ' ;<  a  . scout 
fluid.  Since  balloons  experience  oscillating  rotations  frequently,  it 


appeirs,  at  first,  that  Lamb's  results  can  be  very  useful  in  analyzing 
th'se  rotations.  However,  he  considered  very  small  Reynold's  numbers, 
in  which  case  the  convection  terms  in  the  momentum  equation  of  the  fluid 
can  be  neglected.  Thus,  he  was  able  to  solve  in  closed  form  this  sim¬ 
plified  linear  problem  (the  general  problem,  with  the  convection  terms 
Included,  is  nonlinear).  His  results  for  the  damping  aerodynamic  torque 
and  apparent  .'•‘ditional  moment  of  inertia  depend  linearly  on  the  ampli¬ 
tude  of  osc. i 'at ion.  Notice  than,  in  Equations  12  and  18,  the  dependence 
of  o  and  l_«k  on  u  is  nonlinear,  which  shows  that  the  convection  terms 
are  mcludeo  in  our  treatment.  The  Reynold's  numbers  involved  in  balloon 
rotations  ate  high,  and  Lamb's  results  are  not  applicable. 

Car-i er  and  Di  Prist  (Ref.  7)  considered  the  same  problem  with 
the  convection  terms,  l'ney  formulated  toe  solution  in  the  form  if  an 
Infinite  series  in  ascending  powers  of  the  amplitude  of  oscillation,  and 
they  computed  the  first  linear  term  (a  special  caie  of  Lamb's  solution) 
and  the  next  .(uadratic  terra.:  Their  two-terms  solution  is  accurate  enough 
when  the  following  two  condition*  are  sat(sfied.  The  firat  is  a  condi¬ 
tion  on  the  frequency  of  oscillation  which  can  be  aatiafied  only  by  the 
moat  rapid  balloon  oscillations.  The  second  la,  naturally,  a  restriction 
on  the  amplitude  of  oscilLation  (it  must  be  small)  which  is  not  always 
satisfied  by  balloons.;  Carrier  and  Di  Prima’a  reaulta  are  so  restrictive 
that  thev  can  be  applied  to  the  rotational  oscillation*  of  balloon*  only 
In  special  circumstance* . 
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VXI.  jYKAMICS  of  balloon  rotations 


Now,  that  the  aerodynamic  parameters  of  the  balloon  have  been 
determined,  the  equations  of  motion  may  be  written  for  the  coupled  ro¬ 
tations  of  the  balloon-gondola  system. 

Let  6^  and  62  denote  the  angular  displacements  u f  the  balloon 
anu  gondola,  respectively.  Then  the  apparent  additional  moment  of  in¬ 
ertia  (IA)  and  the  damping  aerodynamic  torque  (D)  of  the  balloon  are 
given  by  (see  Equations  17  and  18):' 


I. 


A 


(20) 


0 


(21) 


A  "dot"  is  used  to  denote  the  derivative  with  respect  to  time.  The  ver¬ 
tical  bars  mean  "alsoiute  value  of".  They  are  necessary,  since  IA  must 
be  positive  and  !)  trust  oppose  the  motion  always.  The  constants  a  and 
fl  are  given  by: 


31  ( pft)l/ 2  a‘ 

(22) 

2.3  (/>yu) 1/2  a4 

(23) 

Let  Ii  ViU  Iy  denote  the  fabric  moments  of  inertia  of  the 
balloon  sad  gondola,  respectively,  and  k  the  effective  torsional  spring 
constant  of  the  suspension  system  between  balloon  and  gondola.  Let  T ^ 
and  ?2  be  tbe  torques  applied  to  the  balloon  and  gondola,  respectively 
The  equations  of  motion  of  the  system  are: 


*1  ®1  4  ®  It  4  £  «lWl/2  *  k  -  ®2>  "  *1  <2*> 

*2  Q2  *  k  “  ^2  V23) 


'XL 


5rtbu-  0  littlr  *nr 


Aerodynamic  effect*  on  the  gondola  have  been  neglected.  For  specified 
applied  torques,  the  angular  displacements  of  the  balloon  and  gondola 
can  be  found  from  the  above  equations. 


A.  SINUSOIDAL  ANALYSIS:  THEORY 

Suppose  that  sinusoidal  torques  (T  l  and  T?)  of  given  ampli¬ 
tude,  frequency,  and  relative  phase  are  applied  to  t.he  sysiem,.  Ulti¬ 
mately,  the  system  will  reach  a  periodic  steady  state.  Due  to  the  non- 
linearities  of  the  system,  the  angular  displacements  (0^  and  $2)  will 
have  an  infinite  spectrum  of  harmonics  in  addition  to  the  fundamental. 
The  spectra  of  0^  and  02  can  be  computed,  but  the  procedure  is  much  mote 
tedious  than  that  of  the  reverse  problem,  that  is  to  say,  for  given 
sinusoidal  0^  and  02  >  compute  the  spectra  of  the  required  T 1  and  T2- 
With  regard  to  the  fundamental  components,  the  two  problems  will  give 
approximately  the  same  resu'ts. 

Let  9}  and  02  be  specified  as  follows. 


91 

=  cos  tat 

(26) 

92 

=  A^  cos  Wt  +  A1 2  sin  (at 

(27) 

where  A[f  A2  and  A 1 2  are  constants..  Then,  it  can  be  shown  that:; 


1a  1-1/2 


(h)Aj)  (l  .11  sin  Ut  4.  0.159  sin  3  ut  *  ....)  (28) 


ojoj1^2  =  -  (WA1)3/2(0.917  sin  ut  -  0.102  sin  3  ut  *.  ....)  (29) 


Substituting  Equations  2b  to  29  in  Equations  24  and  25,  we 
find  that  T2  has  only  the  fundamental  component  ( <0 )  wnile  las  an 
infinite  number  of  componenta  (u  ,  3  cj ,  5  cj  ...).  Let  the  fundamental 
components  nf  t;  and  T2  he  denoted  by,. 


't  ^  »  Tj  coa  h)t  -  T‘ ^  ain  ui 


"X 2  '  c0*  wc  *  T  2  910  Wt 
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where  T^,  T'^,  T2,  and  T'2  are  four  unknown  constants.;  They  are  deter¬ 
mined  from  the  following  four  equations: 


2  3/2  1/2 

Tx  =  (-10  I  +  k)  Al  -  1.11  d  10  '  A1  '  -  k  A2  (32) 


,  .  3/2  .  3/2  , 

I  A  to  A,  -  k  A' 


T'  =  -0.917  yfl  0)  AL 


T2  =  (*w  +  k)  A2  -  k  Aj 


x’2  =  (-  tO~  I2  +  k)  A'2 


Equations  32  to  35  relate  the  fundamental  components  of  the 
applied  torques  and  angular  displacements.  Of  the  seven  quantities  (A^v 
a2>  a'2>  t1i  T'l.  t2>  *n<*  t'2)>  three  must  be  specified  The  remaining 
four  can  then  be  computed.  The  apparent  paradox,  that  we  cannot  specify 
all  the  torques  (Tj_ ,  T'l,  T2,  and  T'2),  is  resolved  when  it  is  noted 
that  the  phase  of  Oj  hat  already  been  specified. 

As  the  harmonic  components  decrease  very  rapidly  (see,  for 
Instance,  Equations  28  and  29),  a  solution,  which  involves  only  the 
fundamental  components,  is  adequate  in  this  highly  nonlinear  problem. 

The  error  is  about  10  percent.; 


B.  alNUSOlDAL  ANALYSIS:  NUMERICAL  EXAMPLE 

As  a  numerical  example,  we  will  consider  the  November  1959 
Project  Strato-lab  Flight  at  cel'ing.  Reasonable  values  for  the  param¬ 
eters  of  the  system  are  as  follows:. 

3  3/2 

Of  *  5.S  x  !0  fr  lbs  sac  ' 


fi  m  4.8  x  103  ft  lbs  sec" 

lx  =  1.5  x  103  slugs  ft2 

1,  .  1,100  slugs  ft2 

k  =  110  fi  lbs 
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We  see  thst,  as  the  frequency  increases,  highs)  torques  are 
required  in  order  to  have  the  same  displacement.  This  is  so,  because 
the  above  frequencies  are  smaller  than  the  antiresonant  frequency  of  the 
system,  which  is  equal  to  about  0.707  times  the  resonant  frequency;  i.e., 
0.22  rad/sec. 


Next  consider  th  >'  case  of  the  same  torque  being  applied  to  the 
gondola.  Then,  the  angular  displacements  of  the  balloon  and  gondola  are 
approximately  in  phase.  They  are  approximately  180  degrees  ou.  of  phase 
with  the  applied  torque.  The  amplitude  of  the  displacement  of  the  bal¬ 
loon  is  again  given  by  Table  VI,  but  the  amplitude  of  the  displacement 
of  the  gondola  is  now  differed-  U'able  VII)  <  We  see  that,  as  the  fre¬ 
quency  increases,  the  displacement  of  the  gondola  becomes  smaller  than 
that  of  the  balloon.  This  is  so,  because,  in  addition  to  the  above  dis¬ 
cussed  antlresonant  frequency  of  the  entire  system,  in  this  case  the 

frequency.  Due  to  the  nonlinearities 
frequency  depends  on  the  amplitude  of 
can  be  shown  that  it  cannot  be  larger 
Obviously,  for  the  larger  frequencies 


gondola  has  another  antiresonant 
of  the  system,  the  value  of  this 
the  applied  torque.  However,  it 
than  (k/I1)l/2  =  0.027  rad/sec. 


of  Table  VII,  the  gondola  is  near  this  antiresonance. 


TABLE  VII 

AMPLITUDE  OF  ANGULAR  DISPLACEMENT  IN  RADIANS  FOR  VARIOUS  VALUES 
OF  AMPLITUDE  (T)  AND  FREQUENCY  (u>)  OF  APPLIED  TORQUE 


sss 

— 

1 

■ 

n 

0.0025 

0.003 

0.01 

o.s 

1 

... 

2 

0.5 

0.9 

1.9 

0.4 

_ 

_ 

1.6 

The  above  angular  displacements  are  of  the  same  order  as  those 
shown  by  the  flight  data.  We  see  that  the  required  torques  *r*  rather 
small.  It  seems  that  they  can  arise  easily  from  aerodynamic  effects 
(the  wake  behind  the  slowly  rising  or  falling  balloon,  etc.).  The  smaller 
of  these  torques  could,  also,  be  caused  by  the  motion  of  matte’’  in  the 
gondola  (telescope,  operator,  navigator,  etc.). 

For  the  itove  frequencies  and  amplitudes  of  the  sppllec  rques , 
the  aerodynamic  damping  on  the  balloon  is  very  small,  but  the  apparent 
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additional  moment  of  inertia  of  the  balloon  is  comparable  to  its  material 
moment  of  inertia. 

The  criteria  for  improving  the  stability  of  the  gondola  depend 
on  the  place  of  application  of  the  torque.  In  general,  the  moment  of 
Inertia  of  the  gondola  will  be,  at  most,  of  the  same  order  as  the  moment 
of  inertia  of  the  balloon.  Then,  when  the  torque  ia  applied  to  the  bal¬ 
loon,  the  stability  of  the  gondola  is  improved  as  k/I,  -*0.  On  the  other 
hand,  when  the  torque  is  applied  to  the  gondola,  the  stability  of  the 
gondola  is  improved  as  k/IL  -♦  oo.  i„  other  wotdi)  for  torque(j  a  d  M 

the  balloon,  the  gondola  must  be  decoupled  from  it;  while,  for  torques 
applied  to  the  gondola,  the  gondola  must  be  well  coupled  to  the  balloon 
These  two  criteria  cannot  be  satisfied  simultaneously.  A  compromise  must 
be  irade,  which  will  depend  on  the  nature  of  the  torques  applied  to  the 
balloon  and  gondola.  Considerable  damping  will,  also,  improve  the  sta¬ 
bility  of  the  system,  but  as  pointed  out  already,  the  damping  in  balloons 
(aerodynamic)  is  very  small. 
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